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ABSTRACT 



In 2000 September we observed a bright X-ray burst from X2127+119 with the Rossi 
X-ray Timing Explorer. This burst has a multi-peaked profile, a peak luminosity of 
cn [ ~6.5xl0^'^ erg s~^, a total integrated energy of ~2xl0^° ergs, and significant photo- 

spheric radius expansion to a maximum extent of i?=118±5 km. From the luminosity- 
^ . temperature relation during the expansion phase we derive estimates for the gravitational 

cn i redshift at the neutron star surface, the corrected Eddington luminosity, and the neutron 

star mass. We detect no slow (~30s) radial oscillations or fast (100-1200 Hz) coherent 
t^ ■ oscillations or QPO during the burst. 

O , The 2000 September event is only the second burst ever observed from this globular 

cluster binary (in M15 = NGC 7078), and it shares many characteristics with the event 
fH , detected by Ginga in 1988 October (Dotani et al. 1990; van Paradijs et al. 1990), the key 

Q^' difference probably being the total amount of material consumed in the thermonuclear 

Q . flash. 

d ■ Subject headings: accretion, accretion disks — binaries: close — stars: individual 



(X2127+119) — stars: neutron — X-rays: bursts — X-rays: stars 
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1. Introduction 

The X-ray source X2127+119 is a bright low-mass 
binary located within 2" of the core of the globu- 
lar cluster M15 (=NGC 7078). It is identified with 
the V--I5 optical star AC211 (Auriere et al. 1984), 
which is far and away the brightest counterpart of any 
globular cluster X-ray source. When folded on the 
17.11-hr optical period (Ilovaisky et al. 1993; Homer 
& Charles 1998), the Hght curve of X2127+119 shows 
a broad symmetrical '^S-hr eclipse at phase 0=0, and 
the largest orbital modulation seen from any X-ray bi- 
nary (1.8 mag in U), but with extremely erratic and 
complex X-ray and optical behavior between phases 
(/)=0.2-0.7. Combining this variability with its low 
Lx/ Lopt ratio of ^20, the natural conclusion is that 
X2127+119 is an accretion disk corona (ADC) source 
observed at a sufficiently high inclination that the 
bulk of the direct X-ray emission is obscured by az- 
imuthal disk structure. However, the detection of a 
very luminous X-ray burst from the source with the 
Ginga satellite (Dotani et al. 1990; van Paradijs et al. 
1990; hereafter "the 1988 burst") not only revealed 
that the compact object must be a neutron star, but 
also that it must be viewed directly with no interven- 
ing obscuration. 

The 1988 burst from X2127+119 had a long dura- 
tion (in excess of 150s), a peak luminosity of 5.1 x 10'^® 
erg s~-^, and a total energy release of G.SxlO'"^ ergs 
(correcting the distance to 10.4 kpc; see below), plac- 
ing it amongst the brightest of the observed Type 
I X-ray bursts. A bright precursor was seen 6 sec- 
onds prior to the main event, and spectral fitting 
revealed strong photospheric radius expansion dur- 
ing the burst, followed by contraction and a cooling 
phase. In the first 30s of this contraction, slow radius 
oscillations occurred (Dotani et al. 1990; van Paradijs 
et al. 1990). 

Per unit mass, low mass X-ray binaries (LMXBs) 
are ~100 times more likely to be found in globular 
clusters than in the Galaxy as a whole, and they are 
almost certainly formed by the tidal capture of com- 
pact objects, an efficient process within globular clus- 
ters because of the high local number density of stars 
(Fabian, Pringle, & Rees 1975; Hut et al. 1992). The 
bursts observed from the 12 known globular cluster 
LMXBs > 10*^^ erg s~^ have a special significance be- 
cause we have relatively accurate information about 
(a) their distance, and (b) the elemental abundances 
of the regions in which they were formed. The M15 



cluster has a low mean metallicity of ^0.01 Solar 
(Sneden et al. 1991; Geisler et al. 1992; Harris 1996), 
although the local metallicity of the accreted material 
in the X2127-I-119 system may, of course, have been 
altered from this by evolutionary effects. 

In this Paper, we report the observation and anal- 
ysis of a second burst from X2127+119, 12 years after 
the first. The burst reported here shares many of the 
characteristics of the earlier event, being highly ener- 
getic with radius expansion, a high peak luminosity, 
and a multi-peaked profile; however, its main energy 
release is of shorter duration, leading to a lower total 
energy, and no precursor is observed. 

2. Observations and Results 

We observed X2127-M19 with RXTE (Bradt, Roth- 
schild, & Swank, 1993) on 2000 August 24 04:36- 
10:31 UT, and then for a series of pointings span- 
ning 2000 September 17 12:02 UT - September 23 
02:36 UT, for an on-source total good time of 150 ksec. 
The data presented here were obtained using the 
"Standard 2" mode with 16s time resolution, and 
an event mode with 16/isec time resolution and 64 
channels covering the full energy range of the PCA. 
The PCA consists of five Xe Proportional Counter 
Units (PCUs) numbered PCUs through 4, with 
a combined total effective area of about 6500 cm^ 
(Jahoda et al. 1996). PCUs 0, 2, and 3 were re- 
liably on throughout the observation, and the data 
shown here are from these detectors. Data were re- 
duced using the RXTE standard analysis software 
included in FTOOLS 5.0.1, and we used XSPEC 
11.0 for our spectral analysis. Where necessary a 
small instrument-related deadtime correction (6%) 
has been applied to the derived PCA fluxes. Errors on 
spectral fit parameters are quoted to 90% confidence. 

The burst was observed during the September se- 
quence, at 2000 September 22 13:47:41 UT, a time 
corresponding to phase 0=O.79±O.O3 in the binary cy- 
cle of X2127-(-119, well away from the partial eclipse. 
In Figure 1 we show the intensity profile of the burst 
in the 2-25 keV range. The profile has a complex 
structure, with an initial fast rise to a first emis- 
sion peak in <0.5s, followed by three further peaks 
in the overall counting rate separated from the first 
by 3.75s, 11.0s, and 21.5s (all such times in this pa- 
per are quoted to ±0.25s accuracy unless otherwise 
stated) . The burst took the source from its persistent 
emission rate of 90 c/s to a maximum of 11,000 c/s. 
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Fig. 1. — The intensity profile of the bright Type I 
burst observed from X2127+119 in 2000 September, 
plotted with a time resolution of 0.1s. 



after whieh it decayed with a I/e-folding time of 30s. 
The effects of the burst are still visible in the X-ray 
spectra 300s after its initial rise. A full timing analy- 
sis was performed on the burst data: no slow (~30s) 
radius oscillations or fast (100 Hz~1200 Hz) coherent 
oscillations or QPO were detected. 

The energy dependence of the burst profile is illus- 
trated in Figure 2. In the 1.9-3.5 keV band the initial 
energy release appears gradual, extending over 6s. At 
the highest energies the burst appears as a bright Is 
pulse, followed by a more gradual 20s rise and an 
almost-symmetrical decay. Such energy-dependent 
double peaking in Type I bursts has long been rec- 
ognized as the signature of photospheric radius ex- 
pansion in bursts reaching the Eddington luminos- 
ity (Hoffman, Cominsky, & Lewin 1980; Lewin, van 
Paradijs, & Taam 1993). 

To study the physical morphology of the burst in 
more detail, we extracted a series of spectral slices 
through the event with a time resolution of 0.25s, and 
performed spectral fits to them over the 2.5-15 keV 
energy range. Our first pass ultimately proved to be 
the most satisfactory: we fit a simple blackbody func- 
tion to the burst spectra, using as background a 300s 
section of persistent-emission data from immediately 
prior to the burst. Experimentation showed that the 
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Fig. 2. — The burst light curve, broken down into 
several bands to show the energy dependence of the 
profile. 



hydrogen column density was poorly constrained by 
the data, leading us to fix it at a nominal value of 
Nh = 7x10^^ cm""^, consistent with previous obser- 
vations (e.g. Sidoli et al. 2000). These fits produced 
acceptable values of x^ throughout the burst. No im- 
provements were gained using more sophisticated fit- 
ting models that include Comptonization of the black- 
body component or modification of the blackbody by 
electron scattering (Nishimura, Mitsuda, & Itoh 1986; 
Nakamura et al. 1989), or a contribution from disk 
refiection of the burst emission (Day & Done 1991). 
We saw no indications of transient narrow absorption 
or emission features in the spectra. 

In Figure 3, we plot the variation of the bolomet- 
ric luminosity, the blackbody temperature fcTf,;,, and 



the derived apparent blackbody radius Rbb through 
the burst as a function of time. (At the temperatures 
and densities found on accreting neutron stars, the 
Compton scattering opacity dominates over the ab- 
sorption opacity, reducing the emissivity below that 
expected for a true blackbody. Thus, this method 
measures the color temperature rather than the effec- 
tive temperature kT^ff. We will return to this point 
in the Discussion section.) In the lowest panel we see 
clear evidence for photospheric radius expansion; the 
apparent radius increases in 5.6s to a maximum en- 
velope radius of i?hf,=118±5 km, followed by an even 
more rapid collapse to i?hfc=29.0±1.5 km over the next 
1.0s. After this point the apparent radius contracts at 
a more steady rate to a minimum value of 8.6±1.0 km, 
20s after the initial burst rise. As is common in radius 
expansion bursts, the measured blackbody tempera- 
ture mirrors the behavior of the derived radius. After 
an initial temperature peak during the first 1.0s of 
the burst, the temperature drops sharply during the 
photospheric expansion phase. During the subsequent 
contraction, the temperature climbs back to its peak 
value of kTbb=2.70±0.13 keV, and then cools expo- 
nentially as the blackbody radius stabilizes. In calcu- 
lating the luminosity we assume that the radiation is 
isotropic and that the source is located at a distance 
of d=10A kpc (Durrel & Harris 1993; Harris 1996). 
The luminosity peaks at a value of 6.5x10^^ erg s~^, 
and over the whole 300s interval where the effects of 
the burst are detectable, the total integrated burst 
energy is 1.85x10^° ergs. 

While it is common practice in burst spectral stud- 
ies to subtract off the continuum emission prior to fit- 
ting, it is well-known that this is a source of potential 
error (e.g. Lewin et al. 1993). If there are two com- 
ponents to the persistent emission, and one compo- 
nent originates on the neutron star or in the boundary 
layer around it, this component of the emission may 
either increase or disappear during the burst itself, 
meaning that the "continuum background" is either 
underestimated or overestimated. This can lead to an 
erroneous decrease in the apparent blackbody radius 
particularly towards the end of the burst, where the 
flux level is low. In addition, for bursts that approach 
or exceed the Eddington limit, the possibility exists 
that the global accretion onto the neutron star itself 
may be suppressed, leaving the burst itself as the only 
source of X-ray emission for the brief duration of the 
burst. 

We investigated both of these possibilities to deter- 




Fig. 3. — The variation of the bolometric luminosity 
(for a distance of 10.4 kpc), blackbody color tempera- 
ture kTbb, and apparent blackbody radius R, derived 
from fits to spectra accumulated through the burst 
with 0.25s resolution. 



mine whether they affect our results. The persistent 
continuum emission prior to the burst (and, in fact, 
through the whole RXTE observation) can be well fit 
using either a straightforward cut-off power law with 
index r=1.63±0.16, Ecut^8.2±1.8 keV, x^=26/36, or 
with a combination of a powerlaw plus blackbody, 
with r=2.42±0.11, fcTb6==1.43±0.16 keV, x^=20/35, 
with the blackbody providing ^^20% of the total flux. 
In the first formulation, there is no contribution from 
the region on or around the neutron star; adopting 
the second, we repeated our series of fits through the 
burst, holding the powerlaw component constant and 
allowing the continuum+burst blackbody component 
to vary. We found that the continuum blackbody 
contribution was a sufficiently small addition that it 
had a negligible effect on Figure 3 and on the physi- 



cal numbers quoted below. (The broadband 0.1-100 
keV spectrum of X2127+119 was found to be well- 
described by a partially-covered powerlaw plus a disk 
blackbody (Sidoli et al. 2000); our spectroscopy re- 
sults are consistent with this model, though the PCA 
lacks the sensitivity to the lowest energies provided 
by the SAX LEGS instrument that would enable us 
to perform a detailed study of partial covering over 
the limited intensity variations found in these data.) 
We next repeated our series of burst fits assuming 
a cessation of persistent emission during the burst, 
and fitting a blackbody to the overall emission with 
only the relevant instrumental PCA background sub- 
tracted. This produced poor fits to the individual 
spectra, with large systematic errors in the residuals, 
indicating that the assumption was invalid and that 
persistent emission does indeed continue unabated 
during the burst. 

3. Discussion 

A comparison of the 1988 and 2000 bursts from 
X2127-I-119 shows that there are strong similarities 
between the two events. Both show multiple peaks in 
their X-ray light curves in count- and flux- space, and 
both show indisputable evidence for radius expansion 
and a qualitatively similar evolution of flux and black- 
body temperature and radius. The bursts have peak 
luminosities of Lmaa;=5.3xl0^^ erg s~^ and 6.5x10^* 
erg s~^ respectively; peak blackbody temperatures of 
kTf,i,=2.87±0m keV and kTi,f,=2.70±0.U keV; and 
developed from comparable levels of continuum emis- 



sion (Fpers=5.0xl0 ^" erg cm 



and 



i^per.=4.2xl0 



^^ erg cm^^ s^^ over the 1-28 keV 
range) (Dotani et al. 1990, van Paradijs et al. 1990, 
and above, assuming d=10.4 kpc). 

The bursts differ in the maximum extent of their 
photospheric expansion - Rbb ^1000 km and Rhb ^120 
km respectively - and in their timescales. The 1988 
burst reached its maximum radius expansion after 
~6s, the end of its sharp radius contraction after 
~18s, and flnal photospheric touchdown (as deflned 
by the simultaneous peak in kThj, and trough in Rbh) 
after ~90s. For the 2000 burst these times are re- 
spectively 5.6s, 6.6s, and ^19s. The e-folding time 
for the decay of the 1988 burst is hard to determine 
from van Paradijs et al. (1990), but they note that 
the burst intensity had decreased to 20% by 160s af- 
ter burst onset. For the 2000 burst, the 20% mark 
is reached after 38s. Largely due to this difference 



in timescales, the total bolometric energy released by 
the bursts differs by a factor of 4: 6.3x10^° ergs for 
the 1988 event, as compared to 1.85x10^° ergs for 
the 2000 event. Although a little weaker than the 
first burst, the burst described here still ranks as an 
unusually energetic event. 

From the fluxes derived earlier we calculate the 
burst characteristics 7 (^Fpers/Fmax) and t {^Eb/Fmc 
where Eh is the total burst fluence or integrated net 
flux across the burst) to be 0.0125 and 125s respec- 
tively for the 1988 burst, and 0.0084 and 28.5s for the 
burst presented above. In their summary of the prop- 
erties of Eddington-limited bursts from a large num- 
ber of sources, van Paradijs, Penninx & Lewin (1988) 
find a strong anti-correlation between log r and log 7 
for Type I bursts showing photospheric radius expan- 
sion. The X2127-f 119 bursts lie on the fine of this 
anti-correlation, but at the extreme end of the (7, t) 
distribution, indicating the intensity and long dura- 
tion of the bursts from X2127+119 relative to those 
from other sources. 

In Figure 4a we present the luminosity-temperature 
(Lboi-kTc) diagram for the 2000 burst, constructed us- 
ing the parameters derived from the fits to the spectra 
out to 100s after the burst rise. The radius expan- 
sion and contraction phases are almost horizontal at 
a luminosity corresponding to the Eddington limit, 
with a slope introduced by general relativistic effects. 
Figure 4a can be directly compared to the equiva- 
lent diagram for the 1988 burst (van Paradijs et al. 
1990, Figure 9); in removing the time domain, this 
representation emphasises the similarity between the 
two bursts. Figure 4b shows the strong inverse cor- 
relation between the apparent blackbody radius and 
temperature; the track defined by the envelope con- 
traction lies almost on top of the expansion track, 
indicating the reversibility of the thermodynamic pro- 
cesses involved. From first principles the luminosity, 
radius, and temperature of the photosphere should be 
related by the expression L = AircrR^T'^. In this case 
the luminosity is not quite constant, and the radius 
and temperature obey the relation R oc T~", with 
a=2.06±0.04 (90% confidence). Both diagrams show 
the radiative cooling curve of the neutron star surface 
after the touchdown of the photosphere. 

While a simple blackbody often provides reason- 
able fits to burst spectra from neutron stars, the en- 
ergy distribution of the 'true' emission is likely to be 
more complex. A systematic difference is expected be- 
tween the blackbody color temperature kTc (—kTbb) 
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Fig. 4. — The data from Figure 3, extended out 
to 100s after the burst rise, and replotted to show 
the (a) luminosity-temperature diagram, (b) radius- 
temperature diagram. Data from the radius expan- 
sion (0s-6s), envelope contraction (6s-19s) and black- 
body cooling curve (19s-100s) are plotted with open 
triangles, filled triangles, and crosses respectively. 
Note the almost-horizontal track in the upper dia- 
gram corresponding to the probable Eddington lumi- 
nosity of the source, and the close correspondence be- 
tween the expansion and contraction profiles in each 
diagram. 



obtained from spectral fitting and the effective tem- 
perature kTeff, due to the dominant contribution 
of electron scattering in the hot neutron star atmo- 
sphere, and the resulting suppression of emissivity 
(van Paradijs 1982). Model atmosphere calculations 
show that the spectral hardening ratio t=Tc/Teff is 
greater than 1, and increases as a function of Tg// and 
of l=Lx/LEdd (London et al. 1984, 1986; Ebisuzaki 
& Nomoto 1986, 1988; Titarchuk 1994; Madej, 1997). 
The neutron star radius observed by a distant ob- 
server, Roo, is related to Rbb by Roo=Rbbt^ ■ Accurate 
knowledge of the behavior of t would place direct con- 
straints on the mass-radius relationship for neutron 
stars. However, there is as yet no convincing agree- 
ment between the various existing theoretical models. 



and in addition the model atmosphere calculations do 
not match the observed variations with luminosity in 
several sources, including X1608— 522 (Penninx et al. 
1989) and X2127+119 itself (see below, and Figure 5). 
If we adopt a representative spectral hardening factor 
of t=1.5±0.1 (e.g. London et al. 1984, 1986; Madej 
1997) for illustrative purposes, our measured photo- 
spheric touchdown radius of i?h6=8.6±1.0 km would 
imply i?oo=19.3±3.4 km. 

The value of the Eddington luminosity measured 
by a distant observer depends upon the mass and ra- 
dius of the emitting object. In the limiting case, the 
Eddington luminosity measured at very large pho- 
tospheric radii, LEdd,oo{R ^ R*)i will be greater 
than the value measured at the point of photospheric 
touchdown, LEdd.oc{R = i?*), by a factor equivalent 
to the gravitational rcdshift at the neutron star sur- 
face, 1 + z^, (Paczynski & Anderson 1986; Damen 
et al. 1990). At intermediate radii, the ratio of 
the Eddington luminosities observed {LEdd,oo{R > 
R*) / L Edd,oo{R — R*)) is a function of both the red- 
shift and photospheric radius, as modified by the 
(possibly varying) spectral hardening (Damen et al. 
1990). Therefore, the photospheric expansion track 
allows us to directly estimate the gravitational red- 
shift factor at the surface of the neutron star, if we 
account correctly for the finite size of the photosphere 
and for the temperature dependence of the scattering 
opacity mentioned above. 

For each data point in the photospheric expansion 
track (Figure 4a), we applied Eqns. 7-9 of Damen 
et al. (1990) to obtain an instantaneous estimate 
of (1 -I- z*), and then took the mean of these de- 
rived values. For this calculation we adopted an 
Eddington luminosity at touchdown of LEdd,oo{R = 
i?*) = (5.08±0.31)xl0'^* erg s~^ and a temperature at 
touchdown of 2.70±0.13 kcV. 

Model atmosphere calculations during the expan- 
sion phase of a burst have not yet been performed, so 
we calculated the redshift for two differing assump- 
tions. First, we assumed that the spectral hardening 
ratio was constant during the expansion, and thus 
that no correction was required for non-Planckian 
shape of the burst spectra. For this case, we ob- 
tained a value of 1 -I- z*=1.17±0.05 (s.d.). In the sec- 
ond derivation we assumed (again following Damen 
et al. 1990) that the dependence of i on the observed 
color temperature Tc^oo is the same during radius ex- 
pansion and contraction as it is during the neutron 
star cooling. We cannot derive absolute values for t 
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Fig. 5. — The variation of the spectral hardening fac- 
tor t^Tc/T(,ff oc Tc/F^-"^^, as a function of the color 
temperature Tc derived from spectral fits. The best- 
fitting straight line for the higher-temperature data 
points is overlayed. 



from observational data alone; we can, however, de- 
termine relative values, and this is all the method 
requires. During cooling, the effective temperature 
Teff, oo is proportional to -F^^^, and thus the harden- 
ing factor t is proportional to Tcoo/F^"^^ (Penninx 
et al. 1989; Damen et al. 1990). In Figure 5 we 
show Tc^oo/F^^ plotted against Tc, oo for the coohng 
track of our X2127-f-119 burst, following photospheric 
touchdown. For temperatures kTc > 1.5 keV, the re- 
lationship can be approximated with a straight-line 
fit. Using this fit to scale the radius estimates that 
enter the redshift calculation (Damen et al. 1990, Eqn. 
9), our derived gravitational redshift factor becomes 
l + z,=1.28±0.06 (s.d.). 

Applying the relativistic correction and solving for 
LEdd and mass (e.g. Lewin et al. 1983, van Paradijs 
et al. 1990), we then obtain, for the constant-t case, 
a corrected Eddington luminosity of LEdd.oo{R S> 
i?*)=(5.94±0.45)xl0'^* erg s~-^ and a neutron star 
mass of Af=2.38±0.18(l + X)Mq, where X is the 
hydrogen content by mass; X=0 for hydrogen-poor 
material and X=0.73 for cosmic composition). For 
the variable-i instance, these values are LEdd,oo{R S> 
i?*)=(6.50±0.50) X 1038 erg s'^ and M=2.60±0.22(l+ 
X)Mq. 

We emphasise that although only statistical errors 
are quoted here, all current derivations of basic neu- 



tron star parameters from burst spectroscopy are af- 
fected by the systematic uncertainties in the interpre- 
tation of the spectral fit temperatures (van Paradijs 
et al. 1990, Lewin et al. 1993). Our own data (Fig- 
ure 5) imply a possible variation in t of ~50% over 
the temperature range 1.5> kTc >2.7 keV which is 
not well reproduced in the model atmosphere calcu- 
lations. The composition of the accreted matter, and 
degree of anisotropy in the burst radiation, are ad- 
ditional key unknowns limiting our ability to obtain 
the mass-radius relationship of the neutron star using 
these methods. 

Nonetheless, while higher than canonical, the val- 
ues of LEdd and M derived here are consistent with 
those obtained by van Paradijs et al. (1990) for the 
1988 burst. It therefore seems likely that the only 
significant difference between the two bursts is in 
the overall timescale of the burst, and thus the to- 
tal amount of matter ignited in the thermonuclear 
flash. The quantity of energy liberated in the combus- 
tion of hydrogen and helium into iron-peak elements 
depends on the ratio of the original mix. Following 
Lewin et al. (1993) and adopting a conversion frac- 
tion of 6x10^8 erg g~^ = O.OOTc^, we calculate that 
the two events respectively consumed l.lxlO^^g and 
2.8xl0'^^g of accumulated material, each equivalent 
to many days of accretion from the secondary. 

We thank Cynthia Hess, Erik Kuulkers, and Tod 
Strohmayer for useful comments and discussions. 
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